INLO-PUB-12/03 



g ■ INSTANTONS AND CONSTITUENT MONOPOLES * 

O , , , 

CS| ■ Falk Bruckmann, Daniel Nogradi and Pierre van Baal 

^ ■ Instituut-Lorcntz for Theoretical Physics, University of Leiden, 

^ \ P.O.Box 9506, NL-2300 RA Leiden, The Netherlands 

Dedicated to the memory of Ian Kogan 

> 

00 , 

, We review how instanton solutions at finite temperature can be seen as 

' boundstates of constituent monopoles, discuss some speculations concern- 

ing their physical relevance and the lattice evidence for their presence in a 
. dynamical context. 

O ■ 

PACS numbers: ll.lO.Wx; 12.38.Lg; 14.80.Hv 

Q^' 1. Introduction 

, Over the last five years there has been a revived interest in studying 

instantons at finite temperature T, so-called calorons ^ The main 
reason is that new explicit solutions could be obtained in the case where 
^ , the Polyakov loop at spatial infinity is non-trivial, necessary to reveal more 

' clearly the constituent nature of these calorons. This asymptotic value of the 

Polyakov loop is called the holonomy. In gauge theories trivial holonomy, 
for which the asymptotic value of the Polyakov loop takes values in the 
center of the gauge group, is typical for the deconfined phase. Therefore, 
caloron solutions with non-trivial holonomy are more expected to play a role 
in the confined phase, still at finite temperature, but where the average of 
the trace of the Polyakov loop is small. In the introduction we start with a 
pedagogical overview discussing monopoles, instantons, and their physical 
significance. The next section discusses the construction of the caloron 
solutions. Sees. ETTI and IT^ are more technical, but tutorial in nature. One 
could skip in particular Secs. l2.3l and l2.41 describing in Sec. Olthe properties 
of the solutions in less technical terms. Lattice evidence for the dynamical 
significance of constituent monopoles is discussed in Sec. [l] 
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1.1. Monopoles 

That monopoles should play a role in describing the constituent nature 
of calorons is in itself not really a surprise, because at finite temperature 
plays in some sense the role of a Higgs field in the adjoint representation. 
However, a gauge transformation, 

<^A^{x) = g{x)A^{x)g-\x) + g{x)d^g-\x), (1) 

shows that Aq does not transform correctly (unless the gauge transformation 
is time independent), due to the inhomogeneous term. Instead, the Polyakov 
loop 



P{t, x) = Pexp ^ Ao{t + s, x)ds^ 



(2) 



transforms as it should, ^P(x) = g{x)P{x)g^^{x). Here (3 = 1/kT is the 
period in the imaginary time direction, under which the gauge field is as- 
sumed to be periodic. We also will consider other gauges, where ^^(x) is 
periodic up to a gauge transformation, in which case the expression for P 
has to be modified accordingly 3 . For example, in the so-called algebraic 
gauge, Ao{x) is transformed to zero at spatial infinity. In this case the gauge 
fields satisfy the boundary condition (Voo to be defined below) 

A^{t + p,x)= VooAf^it, x)V^. (3) 

We will require that the total Euclidean action of these calorons is finite, 
such that the field strength^ 

F^,u{x) = d^A^ix) - d^A^{x) + [A^{x),A^{x)] (4) 

has to go to zero at spatial infinity. It is this that forces the Polyakov loop 
to become constant at spatial infinity. For SU(n) gauge theory this gives 

27rz 

Poo = exp(/?A-)= limP(x), A^ = —U^di^g{^il,^l2....,lln)UQ\ (5) 

>oo p 

independent of the direction and time. Unlike for a Higgs field, however, 
Voo is unitary with determinant 1. Choosing Uq (the constant gauge func- 
tion that brings Voo to its diagonal form) appropriately, the n eigenvalues, 
exp(27ri/Xm), are ordered. The freedom in adding an arbitrary integer to ii„i 
is fixed by ordering the themselves, and requiring their sum to vanish, 

Ya=1 1^-1=0, fJ-l < fJ-2 < ■ ■ ■ < fJ-n < fJ-n+l = 1 + ^1- (6) 



^ Our conventions are Ai_i{x) — eA'^{x)Ta, where e is the couphng constant and = 

-r„, TriTaTb) = -^dab, lTa,Tt] = fatcTc. FoT SU(2) fabc = Sabc and Ta = -i^Ta in 

terms of the famihar PauU matrices. 
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For trivial holonomy Voo is an element of the center of the gauge group, 
Voo = exp(27rig/n)ln with q an integer between and n— 1, hence = Q/n. 

One might now immediately object that a Higgs field that goes to a 
constant at infinity does not have the usual hedgehog form expected for a 
non-Abelian 't Hooft-Polyakov monopole j4j, but this is because the caloron 
solutions are actually such that the total magnetic charge vanishes. The 
force stability of these solutions is based, as for exact multi-monopole so- 
lutions in the Bogomol'ny-Prasad-Sommerfeld (BPS) limit (5j on balancing 
the electromagnetic with the scalar (Higgs) force [HI IZ| • For the caloron the 
difference is in interchanging repulsive and attractive forces. For a single 
caloron with topological charge one, this is because there are n—\ monopoles 
with a unit magnetic charge in the i-th U(l) subgroup, all compensated by 
the n-th monopole of so-called type (1, 1, ... , 1), having a magnetic charge in 
each of these subgroups. This special monopole is also called a Kaluza-Klein 
(KK) monopole 0, as will be explained below. The well-known Harrington- 
Shepard solution has trivial holonomy, and although all eigenvalues of 
Voo are equal, such that there is no spontaneous symmetry breaking, one 
still finds a genuine BPS monopole in a suitable limit (it will be the KK 
monopole that survives). 

In a Higgs theory, switching off the Higgs potential and splitting off a 
square in the energy density, 

- Tr ((A^)' + Sf) = -TV(A$ T Bif T 2Tr (B^A^) , (7) 

exact monopole solutions are constructed using the BPS condition |5j , which 
imposes the covariant derivative of the Higgs field $ to be equal (up to a 
sign) to the magnetic field, A$ = <I>] = ±i?j, where Bi = ^Eij^Fj^. 

One is then left with a total derivative 

TV {BiDi^) = Tr (A(A$)) = 5iTr(Si$), (8) 

whose integral is proportional to the magnetic charge. The integral can also 
be associated to the mapping degree of the map x <I>(rx) (x = a;/|x|) for 
r — > oo. With $ taking values in the algebra, and of fixed length at infinity, 
this gives for SU(2) a map from S'^ to S'^. For a caloron the BPS condition is 
simply a consequence of the self-duality conditions characteristic of instan- 
ton solutions, Ei = DiAo — doAi = ±Bi, with $ = Aq. One might thus be 
tempted to call the constituents dyons, rather than monopoles. In the Higgs 
model the Julia-Zee dyons are constructed by taking proportional to the 
Higgs field $ ^Uj. By a time dependent gauge transformation can then 
be gauged to zero. The resulting electric field is now given by Ei = —doAi 
and is not quantized (and in particular not equal to ±i?j). In pure gauge the- 
ory it makes, however, no sense to separate = DiAq from d^Ai. Gauge 
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invariance requires that they occur in the combination FiQ = D^Aq — doAi. 
The electric field is necessarily fixed and quantized as soon as we interpret 
^0 as the Higgs field. As discussed in Ref. adding a fifth dimension 
for Minkowski time, compactifying the Euclidean time direction to zero size 
(/3 0) allows one to have electric charge as for the Julia-Zee dyon. A 
compactified Euclidean direction is what one also considers in Kaluza-Klein 
theories. It is in this sense that the type (1, 1, ... , 1) monopole is called a 
KK monopole, because it turns out it is static up to a gauge transformation 
that makes one full rotation in the unbroken subgroup when going from 
to (3, thus of lowest non-trivial Kaluza-Klein momentum^. 



1.2. Instantons 

We therefore consider it most appropriate to call the constituents mono- 
poles, as is also clear from Nahm's formalism which provides one of the es- 
sential tools to find these self-dual solutions A3 . The Nahm transformation, 
as well as the Atiyah-Drinfeld-Hitchin-Manin (ADHM) construction [1^ for 
instanton solutions on IR^, form indispensable tools to find the exact caloron 
solutions. It should be added that from the topological point of view Taubes 
showed how to make out of two oppositely charged monopoles a Euclidean 
four dimensional gauge field that has non-zero topological charge J3] . This 
result is more general, since only when minimizing the action in a sector 
with non-trivial topological charge will one find a self-dual instanton solu- 
tion. His construction is based on creating a monopole anti-monopole pair, 
bringing them far apart, rotating one of them over a full rotation (the so- 
called Taubes winding) and finally bringing them together to annihilate, see 
Fig. n The four dimensional configuration constructed this way is topolog- 
ically non-trivial. 




Fig. 1. The topologically non-trivial field configuration is constructed from two 
oppositely charged monopoles, rotating one of them over a full 27r rotation. 



^ Nevertheless, there is a context in pure gauge theories where dyons appear, but for 
this we have to add a term proportional to OFF to the Lagrangian |12J. The electric 
charge is now given by 6/{2n) times the magnetic charge. 
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In the topological charge is related to the winding of a gauge function 
g{x) defined at infinity, which therefore is a mapping from to the gauge 
group. The way g{x) enters is through requiring the Euclidean action, 

S = -^J d'xTr(F2,(x)), (9) 

to remain finite. As before this implies the field strength at infinity to go to 
zero, where A^{x) can be written as a pure gauge, A^{x) = g{x)d^g^^{x). 
For SU{2) a simple parametrization as g{x) = a^{x)a^ in terms of a unit 
vector and unit quaternions ci^ = (l2,iT) (ct^ = makes this winding 
most transparent as the mapping degree of a{x) at infinity. In the simplest 
case of degree 1, a{x) = x/\x\ is a 4-dimensional hedgehog. The relation to 
Taubes' construction is using the fact that can be viewed as a twisted 
product of (the Taubes winding) and 5^ (the hedgehog formed by the 
Higgs field), the so-called Hopf fibration PH [TH] . 

In the case of periodicity in the imaginary time direction as it occurs 
for calorons, let us transform ^0(2;) — > everywhere. This can be done by 
a time dependent gauge transformation U{x), which in general will not be 
periodic. Actually P{t; x) itself is the gauge transformation that relates the 
gauge field at t + /? to that at t in the Aq = Q gauge. Since P{x) goes to a 
constant (T^oo) at spatial infinity, this provides a non-trivial mapping from 

(as IR^ compactified at infinity) to the gauge group. The topological 
charge is precisely its winding number. 

As in the case of the monopole energy density, we can rewrite the action 
density in terms of a square and a boundary term, 

TV(i^2j ^ iTr(F^,±F^,)2^TV(F^,F^,), T,{F^,F^,) = d^K^, 

K^, = 2e^^^pTv{A^da,Ap + lA^A^Ap), (10) 

with F^y = ^e^iyapPa/s the dual field strength interchanging electric and 
magnetic components. Hence, for self-dual solutions at finite temperature 

S = ^J {Ko{0,x)-Ko{P,x))d^x = ^J TT{P{x)dP-\x)f, (11) 

using that in the Aq = gauge A{j3,x) = ^A{0,x), expressing the result 
in a compact differential form notation. Note that the integral is invariant 
under smah deformations of P, using 6{PdP-^f = 3dTy:{PSP-^{PdP-^)^), 
and therefore should be proportional to the winding number of P. A map 
of degree k can be obtained by multiplying k maps of degree 1 (the inverse 
gives a negative mapping degree). Each of these maps of degree 1 can be 
deformed at will, as the integral is invariant under continuous deformations, 
and in particular can be arranged to be the identity except for a small 
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region. Choosing these regions to have no overlap, the integral is easily seen 
to be proportional to k. To fix for SU(2) the constant of proportionality 
we may take U{x) = ((1 — |xp)cro + 2x • c?)/(l + for the map of 

degree 1 related to stereographic projection from to IR^. For SU(n) 
we first deform the map to lie in an SU(2) subgroup. This then gives 
the celebrated result that self-dual solutions with given topological charge, 
k = {16tt'^)~^ J d^xTr (^F^^{x)F^^{x)^ , have an action equal 87r^|A:|/e^. 

This is also a convenient setting to understand why in the limit of zero 
temperature, /3 — > oo, an instanton corresponds to vacuum to vacuum tun- 
neling. Finite action requires the field strength to go to zero at \t\ oo, 
but at the same time, the field at t ^ oo is related to the field at t — > — oo 
by a topologically non-trivial gauge transformation. Strictly speaking, we 
need to identify gauge field configuration that are gauge equivalent. Instead 
of having multiple vacua we can alternatively say there is one vacuum, with 
the field space being non-contractible. This is quite analogous to consider- 
ing a periodic quantum mechanics problem, which could be reinterpreted as 
quantum mechanics on a circle. Tunneling now corresponds to penetration 
of the wave function in the classically forbidden region, to reach back to 
the vacuum going around the circle in either direction. Essential is that the 
support of the wave function, the region where it is non-vanishing, becomes 
sensitive to the non-trivial topology of the configuration space |3j . 




Fig. 2. Shown is a typical example for the action density (on equal logarithmic scales 
at the time where it is maximal) of an SU(2) caloron with non-trivial holonomy 
^2 = = 0.125 (for, from left to right, p//3 = 0.8, 1.2, 1.6). 



The tunneling path in non-Abelian gauge theories is of course described 
by the one parameter (t) family of gauge fields Ai{t, x) (in the = gauge), 
and a finite action means we have to cross a potential barrier when going 
around the non-contractible loop (or going from vacuum to vacuum). From 
Taubes' argument it is already clear that this intermediate configuration can 
be associated to a monopole-antimonopole boundstate, which is made more 
precise by the caloron solutions to be described. At zero-temperature this 



7 



is less obvious, since there monopoles form close pairs, cmp. Fig. |21 (left). 
They behave very much like virtual particles, that can only be created for 
a short period of time, thereby explaining the instantaneous character. At 
finite temperature, the monopoles can be separated much more easily due 
to the interactions with the thermal bath. At high temperature, however, 
they will be suppressed due to the Boltzmann weight, and so monopoles 
(like calorons) are expected to be dilute. At low temperature instantons 
form a more dense ensemble, possibly leading to monopoles to be dense as 
well, particularly when instantons overlap. One may then have some hope 
that the confining electric phase could be characterized by a deconfining 
magnetic phase, where the dual deconfinement is due to the large monopole 
density, in a similar spirit to high density induced quark deconfinement. It 
would offer a possible alternative to existing scenarios. 

Instantons describe virtual processes and are quite often discussed in the 
context of the semiclassical approximation. For a one-dimensional particle 
with mass m in a positive potential V we may again split off a square, 

L{t) = ^mx^{t) + V{x{t)) = \m (x{t) - ^1V{x{t))ln^ + x{t)^2mV{x{t)), 
Sci = j dt ±{t)^2mV{x{t)) = j dx\l2mV(x), (12) 

where we typically integrate between the classical turning points, related 
to the WKB expression for the wave function in the classically forbidden 
region, exp(— /^^ dy\/2mV{y)/ti). In a double well this leads to a tunnel 
splitting proportional to exp(—Sci/h). One problem in the theory of strong 
interactions is that the effective coupling tends to become too big for large 
instantons. Instantons can have an arbitrary size p, due to the classical 
scale invariance of the theory. This classical scale invariance gets broken by 
the regularization and one is left with a scale dependent running coupling 
constant after renormalization. This causes a problem when integrating 
over the scale parameter in the one-instanton tunneling amplitude given 
by the celebrated result of 't Hooft [17 , dW oc dpdx'^p'^ exp (—Svr^/e^ (/?)). 
Actually, for calorons with non-trivial holonomy and p > (3, it is more 
natural to associate p with the distance between constituent monopoles 
(irp^/P for SU(2)), and this may help alleviate the problem one encounters 
when dealing with large scale instantons. This is also one way to understand 
why there will still be calorons with p arbitrary, despite the fact that P fixes 
the scale^. At zero temperature, a large p leads to a low energy barrier along 
the tunneling path, and at some point this will no longer describe a virtual 



^ Remarkably, it can be shown |18l I19| under some mild conditions, that any four 
dimensional manifold will have instanton solutions with arbitrary p. 
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process and the semiclassical approximation will break down. Nevertheless, 
the instanton liquid model has been very successful in describing much of 
the low energy phenomenology, in particular for chiral dynamics and aspects 
related to breaking the axial U(l) symmetry. We refer to the reviews by 
Schafer and Shuryak [20] and by Diakonov 21 for more details. 

1.3. Fermion zero-modes 

In the instanton liquid model the interaction of instantons with fermions 
plays an important role. This is because instantons have a remarkable in- 
fluence on the spectrum of Dirac fermions, it namely implies the presence 
of zero eigenvalue solutions with fixed chirality (the so-called chiral zero- 
modes). These chiral zero-modes will play an important role in the Nahm 
transformation and in studying the properties of the calorons with non- 
trivial holonomy, so we wish to mention some of the interesting and far 
reaching physical consequences. This is most easily discussed in terms of 
the spectral flow of the Dirac Hamiltonian along the tunneling path. Its 
spectrum is of course gauge invariant, and this implies that the energy lev- 
els at t = and t = (3 are identical. The gauge field provides a smooth 
interpolation between these, which leads to each of the energies to be a 
continuous function of t, the so-called spectral flow, see Fig|31 If we draw 




Fig. 3. Schematic representation of the spectral flow. 

the energy level at zero, the number of crossings is clearly a conserved 
quantity, it cannot change under continuous deformations of the gauge field 
background (the background need not be a solution of the equations of mo- 
tion). It is also clear, by putting instantons in a row, that the number of 
crossings is proportional to the topological charge, and not surprisingly it 
is actually equal to it. The necessity of such a crossing is related to the 
existence of a chiral zero-mode, whose number for fermions in the funda- 
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mental representation is equal to the topological charge, as follows from the 
Atiyah-Singer index theorem . The argument is most simple in the case 
of zero-temperature, where the zero-mode at t — > ±00 has to behave as 
exp(—tE±), where E± is one of (free) fermion energies. Clearly this would 
only give a normalizable zero-mode if > and E- < 0, which forces a 
crossing! At this point the Dirac vacuum is degenerate and particles or anti- 
particles are created with definite chirality. Compatible with the breaking 
of the axial U(l) symmetry, the divergences of its current is proportional 
to the topological density Ffj,i,Ffj_y. It gives rise to the so-called 't Hooft 
interaction 17 , which together with a finite density of fermion modes at 
zero eigenvalue (required for the Banks-Casher mechanism ^3] to work) 
play an important role in the spontaneous breaking of chiral symmetry (or 
soft breaking in case of small up and down quark masses) j2Uj . Finally, 
the spectral flow also makes it easy to understand the origin of the axial 
anomaly, which occurs due to the need to regularize the theory. If cutting 
off modes with E < —A, where A is the ultraviolet cutoff, we see that the 
spectral flow leads to the fact that modes we had removed in the trivial 
vacuum reappear due to the spectral flow at the vacuum equivalent to it 
by a topologically non-trivial gauge transformation. That the violation in 
conserving the axial U(l) charge is proportional to the topological charge 
is in this setting simply a consequence of the number of crossings in the 
spectral flow. This is the celebrated infrared-ultraviolet connection, and 
also makes it understandable why the anomaly is robust (fully determined 
by the lowest order result in perturbation theory). 

2. Construction of solutions 

Let us start with the well know SU(2) Harrington-Shepard solution ^ 
for the caloron with trivial holonomy. In that case one simply takes a 
periodic array of instantons parallel in group space (Poo = 12)1 placed at 
a(p) = (ao + pf3, a) for integer p. The solution is found in terms of the 
't Hooft ansatz |24 , which in its general form is given by 

Af,{x) = ^fif,^d^ log (t){x). (13) 

Here r/^,, = ^{af^a^ - a^a^) = ir/'^^Ta and fj^^ = \{a^ay - Oyfj^) = ifi^^^Ta, 
are the self-dual and anti-selfdual 't Hooft tensors. Substituting this in the 
self-duality equation, F^y = F^^, one finds this is a solutions if and only if 
5^(/>(x) = 0. Therefore (p = 1 + J2i=i pfl^ ~ where are the four 

dimensional locations and pi the sizes of the k instantons'*. A singularity 

Using conformal transformations a generalization of this ansatz including non-trivial 
color orientations exists |25| . but it will only give all possible solutions for charge 2. 
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at X = can be removed by a gauge transformation. For the Harrington- 
Shepard solution, taking = a(p) and pp = p one can perform the sum 
over p £ Z to find P , 

. / N , 7rp2sinh(27rr//3)/(/3r) _ ^, 

cosh(27rr/p) — cos(27r(t — ao)/p) 

The case of non-trivial holonomy cannot be treated in the same way, 
because we will have to sum over a periodic array of instantons that has a 
color rotation by Voo when shifting time over p. This can only be treated 
within the full ADHM ansatz. The rest of this section is more technical and 
could be skipped, although the short tutorials on the ADHM construction 
and the Nahm transformation in Sees. I?!T] and are still recommended. 



2.1. ADHM formalism 

The SU{n) ADHM construction for charge k instantons 14_ starts with 
a k dimensional vector A = (Ai, . . . , A^), where A| is a two-component spinor 
in the n representation of SU (n) (i.e. A is a n x 2/c complex matrix) and a 
2k X 2k complex matrix B = a^® B^ (each is a hermitian kxk matrix). 
These are combined to form a {n + 2k) x 2k dimensional matrix A(x), 
which has n normalized eigenvectors with vanishing eigenvalue, combined 
in a complex matrix v{x) of size (n -|- 2k) x n, 

^(^)=(i?t))' ^(^) = ^-^l'fc' At(x)t;(x) = 0. (15) 

Here the quaternion x = x^o"^ (a 2 x 2 matrix with spinor indices) denotes 
the position and v{x) can be solved explicitly in terms of the ADHM data, 

<^) = {~1^))<I''^'^^ u{x) = {B^-x^)-^X\ 0(x) = l„+nt(x)n(x). (16) 

The square root cj?-/'^{x) is well-defined because 4){x) is a positive n x n 
hermitian matrix. The gauge field is now given by 

A{x) = A^{x)dx^ = v\x)dv{x) (17) 

and to check its self-duality it is best to use form notation. For the field 
strength two form F{x) = ^F^^{x)dXfj_ f\dxy = dA{x) + A{x) f\A{x) we find 

F = d{vUv) + vUv A vUv = dv'^ A {1 - V (g) v^)dv. (18) 

We note that v ® v'^ projects to the kernel of A^^. Assuming that A has 
no zero eigenvalues such that A^A is invertible, we find that 1 — f f ^ = 
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A(A^A)^^A^^. Indeed, when acting on elements in the kernel of A^, left- 
and right-hand side are equal. Any vector in the orthogonal complement of 
this kernel can be written as Aw, since < v,Aw >=< A^v,w >= 0, such 
that left- and right-hand side are again equal. Hence 

F = dv^A A {A^A)-^AUv = dv^bdx A (AtA)-^^^^^^^^;, (19) 

where we use the fact that A^dv = {dA^)v = dx^b^v, with = (0 I2 ^ Ik) 
as a 2A; X {n + 2k) matrix^. It might seem this does not help that much, 
but remarkably, using dx = dx^a^ we find dx A dx'^ = "q^i^dx^ A dxy, and 
since rj^y is self-dual, we are done. Not quite so yet! We had to take dx ® Ik 
through (A''"A)~^ and this is only possible if 

(At(x)A(x))-i = l2®/., (20) 

where fx is a hermitian kxk matrix. This condition, stating that A^ (x) A(x) 
is invertible and commutes with the quaternions, is what is known as the 
quadratic ADHM constraint. It has reduced solving a set of non-linear 
partial differential equations to solving a quadratic matrix equation. 

To construct a charge k caloron with non-trivial holonomy [201 5 place 
k instantons in the time interval [0, performing a color rotation with Voo 
for each shift of t over cmp. Eq. (PJ- This is implemented by requiring 
(suppressing color and spinor indices and scaling /3 to 1 in this section) 

\k+k+a = "PooApfe+a, Bpk+a,qk+b = Bpk-k+a,qk-k+b + O'O^pq^ab- (21) 

Solutions to these equations are parametrized by (a and Ap^, 

\k+a = 'P^Ca, Bpk+a,qk+b = PO'oSpqSab + ^p-g, (22) 

with A to be determined by the quadratic ADHM constraint. 

2.2. Fourier-Nahm transformation 

We introduce the n projectors Pm on the mth eigenvalue of Vcxn, such 
that Voo = E^e^-^/^-P^ and Xpk+a = Y^m^^'^""^^'^ PmCa- Fourier transfor- 
mation now leads to 

27TiY,^^'P'Af = A-\z), 5]e-2-^'^Apfc+, = 5]5(z-/x^)P^Ca. (23) 

p p m 

^ The notation dx^'b'^ niay be a bit misleading, it should be read as {dx^ (g) lfc)6^. 
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Here is again a two-component spinor in the n representation of SU{n) 
and A°'^{z) = a^A'^{z), with A^{z) an anti-hermitian k x k matrix. In 
terms of the latter 

E B,,^a,,k^,{x)e'-^^P^-'^^'^ = ^-^^^Dt(z'), (24) 

dz 

which is the Weyl operator (positive chirahty Dirac operator) for the V{k) 
gauge field ^^(z)— 27rix^lfc defined on the circle, zG [0, 1], i.e. with period- 
icity 1 in case (3^1). The quadratic ADHM constraint now reads 

ai Ifc, Dliz)D,iz) + 47r2 ^ <5(z - ;U^)Cl^'mCfe] = 0. (25) 

m 

Introducing 

27rClPmCb^t2S^-r-pZ\ (26) 
this leads precisely to the so-called Nahm equation (13j . 

^Aj{z) + [Ao{z),Aj{z)] + ^ejM[Mz),Mz)] = 2vriE'^(^-Mm)Pil. (27) 

Note that the left-hand side is the difference between the magnetic and 
electric field for Af^{z), with the right-hand side a violation of self-duality. 

Although we do not want to go into too much detail here, it is instructive 
to discuss the standard setting of the Nahm transformation |13[ I27j . One 
starts from an SU(n), charge k self-dual gauge field on IR^ with periods in all 
four directions, some of which may be infinite or zero (in the latter case effec- 
tively leading to a reduced dimension) . This extends to a family of self-dual 
U(n) gauge fields^ when adding —2-Kiz^ln to the SU(n) gauge field A^(x). 
The index theorem now guarantees there is a family of k chiral zero-modes 
for the Weyl equation D\^{x]z) = —a^D^{x;z)'^{x]z) = 0, and this al- 
lows one to construct the gauge field A'^{z) = J d^x^'"(x; z)^9^'^(x; z)/dz^. 
This dual U(A;) gauge field can be shown to be self-dual with charge n, and 
is defined again on IR^, but with its periods inverted. Remarkably, one can 
then perform the transformation again, and come back to the original gauge 
field 1271. 

The ADHM construction performs precisely this last step. For instan- 
tons on IR^ all periods are infinite and the Nahm transformation reduces 
self-duality to algebraic equations. Singularities may appear (except when 

® One easily checks it does not change the field strength Ffi„{x). 
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all periods are finite), as we have seen from our analysis using Fourier trans- 
formation of the ADHM data and like in IR^ are related to introducing A, 
which encodes the asymptotic behavior of the zero-modes. 

One final result is crucial to appreciate the strength of this formalism, 

dId, = -I2 ® Dl{x- z) - f]^^ [Df,{x; z),D^{x; z)], (28) 

which uses the fact that a^ai, = 6^u'^2 + ^^lu- Since [Df^{x; z), Di^{x; z)] = 
-F^i/(x) and the contraction of an anti-self dual tensor (r/^^i/) with a self-dual 
tensor (-F/^jy) vanishes, self-duality implies that D\Dz = — 12 ® D^(x;z), 
which therefore commutes with the quarternions. This is completely anal- 
ogous to our discussion for the ADHM construction, and proves that the 
Nahm gauge field A^{z) is self-dual. Doing the Nahm transformation for 
the second time leads us to perform the same calculation as in Eq. this 
time for the dual Weyl operator D^, and we see that the quadratic ADHM 
constraint is fully equivalent with stating that the dual gauge field is self- 
dual. The only complication is the possible presence of singularities when 
some of the periods are infinite, in particular for the ADHM construction 
somewhat hiding this profound relationship to the Nahm transformation. 

2. 3. Some explicit formulae 

One might think all of this is just pushing the same problem around, 
but as we noted before, there is a dramatic simplification due to the di- 
mensional reduction. In addition, for topological charge 1 the dual gauge 
field is Abelian, further simplifying the Nahm equation. In that case all 
the commutator terms in Eq. H27() vanish, and A^{z) is piecewise constant, 
only jumping at the singularities. It is this that has allowed us to make 
progress in finding explicit solutions, together with the "magic" formulae 
for gauge field, action density, fermion zero-modes and density found previ- 
ously within the ADHM formalism pglOU] . 

A,{x) = icp'/\x)f]l,d,^,{x)^y\x) + ^,[ct>-'/\x),d^ct>'/\x)], (29) 

^^,(x) = (2vr)-i {<P'/Hx)Xd,f,a^e)[^, ¥,jixr^^ix) = -i27rr'dlfj^, 

where I = 1, . . . ,k labels the zero- modes (i, / are the gauge and spin index), 

0(x) = {In - 0o)"^ , = A(fT^®/^)At, e = a2 = ir^. (30) 

The gauge field for SU(2) further simplifies to Afj_{x) = \4){x)fi^^j^dy(j)j{x) 
(which could be viewed as a generalized 't Hooft ansatz), since in that case 
(t>{x) as a 2 X 2 matrix is a multiple of I2. For the action density one 
finds [ini, 

TiFl, {x) = -dldl log det (31) 
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It is thus very convenient to first find the matrix fx, as defined in Eq. (|2()|) . 
For calorons, after Fourier transformation, it is replaced by the Green's 
function fx{z,z'), which satisfies the equation j26j 



D^{z;x) 



^ +^5]<5(z-//™)5^|/,(z,z') = M(z-/)- (32) 



2TTi 

The gauge field, as determined through 0^(2;), see Eq. (|29|) . is read off from 

</.^(x) = P^nCa^Jfif^m,l^m')ClPra', (33) 

whereas the fermion zero-modes satisfying a generalized boundary condition 
#^(t + l,x) = e2'^*^7^oo^"(i,^), cmp. Eq. ©, are given by PMI^ 

l'^(x) = (2^)-V'/'(^)E^-Cfe^M^^^/f(Mm,^)- (34) 

m 

For 2 = 5 this gives the usual anti-periodic boundary conditions for fermions 
at finite temperature, but the general z dependence will turn out to be 
extremely useful as a diagnostic tool. Finally, the zero-mode density reads 

^'»(x)t^^(x) = -{27r)-'dlff{z, z). (35) 

To compute fx{z,z') we note that we can always transform to a gauge 
where Aq{z) = 27ri^o is constant, after which g{z) = exp(27ri(4o — xolk)z) 
transforms Aq{z) — 27rixo to zero. This turns Eq. (|32() into 

i^-^ + V{z;x)^fx{z,z') = 4TrHk6{z-z'), (36) 

with fx{z,z') and V{z;x) given by 
fx{z, z')=g{z)fx{z, z')g\z'), V{z; x) = An^ {z; x)-F27r^ 5{z-fim)Sm, 

m 

Rj{z; x) = Xjlk - {27riy'^g{z)Aj{z)g\z), Sm = 9{fJ'm)Smg\fJ-m)- (37) 

Periodicity is now only up to the gauge transformation ^(1). In particu- 
lar when R'^{z;x) is piecewise constant explicitly computing the Green's 
function becomes doable. Nevertheless, in general terms one finds 

fx{z, z') = 47r2 [W{z, zo) (^e{z' - z)l2k - {l2k - J^z,)-^) W-\z', Zo)}^^ , 
W{z2,zi)^PexpJ^'' (y(z;x) o') -^-o = 5^(1)^^(^0 + 1, ^o)- (38) 
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where the (1, 2) component on the right-hand side for fx{z, z') is with respect 
to the 2x2 block matrix structure. This equation for fx{z,z') is vahd for 
z' € [zjZ + 1], but can be extended with the appropriate periodicity. We 
can now also find an explicit result for the action density 

TtF^^,{x) = dldl log^(x), ^l,{x) = det (ie--"'"(l2fc - ^.o)/^) • (39) 
Note that zq can be chosen at will, e.g. for zq = fim + we find 

•^Mm = Tmlim-l ■ ■ ■ T2HiTig^l)HnTnHn-l ■ ■ ■ Tm+lHm, (40) 

The main pay off to have these expressions in terms of possibly unknown 
solutions to the Nahm equations, is that it allows us to look at the far 
field limit. When constituent monopoles are well separated, the charged 
components of the field in the core decay exponentially, being left with the 
Abelian fields in the far field region. This is related to the high temperature 
limit, in which case the cores collapse to zero size. 



2.4- Limits and special cases 
To extract exponential factors it turned out to be convenient to define 



/m(^)= Pexp 



±27r f Rt{z)dz 



Rt{zf±^^^Rt,{z) = R'[z;x). (41) 



Since R{^z\ x) xt^ for |x| — > oo, we find in this limit that R^{z) — > l^llfc 
and f^iz) exp (±27r|x|(2; - Hrn)tk)- For z,z' G [fim, fJ-m+i] we can now 
write W{z,z') = Wm{z)W-^{z'), see Eq. (jSHI), with 

'^^^^^'> - [27rR+{z)f+{z) -27rR-{z)f^{z) ) ■ ^^'^ 

This has allowed us to show that (ff stands for far field limit) 

f^{fim,IJ"m) = 2vri;„\ and f^{z,z) = 7ri?„^(z) for fj,rn<z<fJ.m+i, (43) 

where = -Rm(^m) + Rm-i(t^rn) + Sm and Rm{z) = |(-Rm(-2) + Rmiz))- 
For this to hold x has to be far removed from any constituent. In addition 
mj fJ-m+i) — 0, which implies that for SU(2) and SU(3) only the Abelian 
components of the gauge field survive in this limit. 

Useful is also the so-called zero-mode limit (zm), which assumes x to be 
far removed from any constituent monopole other than of type m, distin- 
guished by their magnetic charge and mass Sn'^i'm, with = fJ-m+i — fJ-m, 
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see the next section. In this hmit one finds up to exponential corrections [M\ 
for IJ-m< z' < z < iim+i {fx{z\ z) = f^{z, z') for z' > z) 

fT{z,z') = vr [f-{z)f-{^^^+l)-^- x (44) 

with Z- = tk- 2E~^i?m-i(/"m) and Z+ = Ifc - 2S;;^iim(/Um). One con- 
cludes that the zero-modes (see Eq. (|5H) and Eq. (|SS1)) are locahzed expo- 
nentially to the mth constituent monopoles (for z away from the interval 
boundaries). 
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Fig. 4. Constituent locations y based on Eq. (|45|) as a function of ^ = = —^2 for 
y2 — —yi = 1^1 = = 5 and pi = p2 ~ 2. Dashed versus full lines distinguish the 
magnetic charge of the constituents. The dotted lines apply to Figs. 13 

Explicit solutions are found when R{z;x) = x — Y{z) is piecewise con- 
stant, i.e. Y{z) = Yrn for z G [fim, ^J-m+i], in which case |^| = = Rm- 
Apart from k = 1, this is so for a class of axially symmetric multi-caloron 
solutions 26 , in which case 3A: eigenvalues of Ym give the locations of the 
type m constituent monopoles. The expression for the zero-mode in that 
case can be used to show that in the high temperature limit the zero-mode 
densities reduce to delta functions localized at these constituent locations, 
thus establishing that the constituents in this limit are point-like monopoles. 
Particularly simple axially symmetric solutions are found for SU(2), 

CiI = PaSiI, Y^ = e{6abia + yraPaPh), ^^12/1 + ^^2^2 = 0, > 0. (45) 

For charge 2, choosing pi = p2 = 2 and equal mass constituents (z^i = 1/2 = 
I) the constituent locations y are plotted as a function of ^ = ^1 = —^2 in 
Fig.|lJ We note that opposite charges are found to alternate. The properties 
of these solutions will be further discussed in the next section. 
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More generally, for SU(2) and charge 2 we may borrow from the study of 
monopoles ^1 EH the general solution of the Nahm equations for z G [/ii , 

Aj{z) = 2^ig\z)h) (ajt2 + Vnjbfb{4TTVz)n)hg{z), (46) 

where fj{z) are the Jacobi elliptic functions 

cnk(^;) cnk(-z) cnk^[z) 

Here a is the center of mass for monopoles of a given type, TZ and h are 
arbitrary spatial and gauge rotations, I? is a scale parameter, and < k < 1 
(k' = Vl — k^) playing the role of a shape parameter, as will be discussed 
in the next section. In general all these parameters differ on the second 
interval where in addition z is shifted to z — ^, but they are to be related 
through the discontinuities in the Nahm equation, Eq. (|27j) . This tends to 
be rather restrictive, and is the point where the construction for calorons 
deviates from that for static monopoles. 

Quite remarkably, although R{z;x) is no longer piecewise constant, the 
function Tri?~^(z) is independent of z £ [nm, fJ'm+i]- This is a highly non- 
trivial consequence of the Nahm equations (from the point of view of in- 
tegrability, it gives a constant of motion). The physical significance here 
is that it is directly related to the zero-mode density (summed over the 
zero-modes) in the high temperature limit, see Eq. H35|) and Eq. (|43|) . 

= -dfVn^ix), V^(x) ^ (47r)-iTrii„i(z). (48) 

a 

Miraculously we have been able to calculate Vmix) exactly jSJ, from which 
we will be able to draw conclusions on the pointlike nature of the con- 
stituents. 



3. Properties of solutions 

We start our tour of SU(n) caloron solutions with those of charge 1. In 
this case the action density has a particularly simple form [^H] 

-iTVF^i(x) = -15292 log V;(x), ^{x) = itr(^„ • • • ^1) - cos(2^V/3), 

^ =_}_[rm \pm+i\\ fcosh(2iTiymrm/(3) sinh(27rj/mrm//3) \ .^gx 
~ -Tm \ / V sinh(27rf„r„/^) cosh(27rfmrm//3)/ ' 

Here = iJ.m+i-IJ'm, Tm = \rm\, Tm = X - Vm and pm = Vm-ym-i-, where 
ym {Vn+m ^ Vm) are the constituent locations''. A few typical examples for 

^ This can be derived from Eqs. 1391401 using Sm ~ Sm = \p?n\ (see Eq. 12611 ) and 

= B,„+ir„+ii?,„B~\ where = B?^ = 0, B^ = 1 and B?^ = 4Tvr^. 
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SU(2), where |pi| = vrp^//? in terms of the instanton scale parameter /O, are 
shown in Fig. Inland for SU(3) in Fig. [S] (we apologize for a somewhat awk- 
ward choice of conventions in naming the vectors pm and scale parameters 
Pa = iCaD- From this it is already clear that there are n lumps, centered at 
ijm and that when well-separated they are static and spherically symmetric. 
The self-duality then guarantees each lump is a basic BPS monopole, which 
can be seen to contribute Sn'^i'm to the action, correctly summing to a total 
action Stt^. Interesting properties on the geometry of the moduli space and 
alternative approaches can be found in Refs. 051 OS] , 

It is instructive to provide further evidence for the monopole nature of 
these constituent lumps. For this we take the far field limit (ff) of Eq. (|49|) . 
In this limit x is assumed to be far from all constituent locations ym and 
we find 

m=l ^'^"^ 

Note that there is no longer any time dependence and the far field limit is 
equivalent to the high temperature limit (3^0, where the monopole mass 
becomes infinite and the non-Abelian core collapses to zero size (the range 
of the exponentially decaying charged fields shrinks to zero). As we will 
argue, using A^{x) as the gauge field for the basic self-dual Abelian Dirac 
monopole, one finds in a suitable gauge for the Abelian far field embedded 
in SU(n), Afix) = 5abA'l,{x), with i-(f) = i^°(f™) - ^A]i{rm-i), where 

E^{x) = BP{x) = are electric and magnetic fields of this self-dual 

Dirac monopole. One verifies*^ that indeed —\TiF^^{x) = —^Ofdj log V'^(x) 
when substituting this gauge field. 



Fig. 5. Action densities for the SU (3) caloron with (/ii, /i2, /^a) = (^17, —2, 19)/60 
and (i^i, 1^2, 1^3) = (0.25, 0.35, 0.4) at t = in the plane of the three constituents for 
1//3 — 1.5, 3 and 4 (from left to right) on equal logarithmic scales. 



-3 -* —3 -* \ 

ri — r„ r2] 



^ It is best to first check —^dld^loglr-^ ""^rj ^{ri + r2 + \p2\Y] = (r^ 
relevant for SU(2). The exponential terms in Eq. 150^ do not contribute, except 
for J(r,„)87r^!/m//3. However, in the core Eq. 1501 should not be used since the 
singularity in ^^^{x) is simply due to approximations made outside the core. 
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To give a little more insight in why the gauge field takes the above form, 
we remark that another way to define the location of an SU(2) monopole 
is to find the zeros of the Higgs field, or for SU(n) to find where two of 
its eigenvalues coincide. At these points the broken gauge symmetry is 
partially restored to include an unbroken SU(2) subgroup. In the case of 
the caloron we of course need to find coinciding eigenvalues of the Polyakov 
loop^, P{x) |S7j. The gauge field A"^{x) defined above is associated to 
the Abelian component of the embedded SU(2) monopole associated to the 
unbroken subgroup. The Polyakov loop therefore is a useful diagnostic tool 
particularly when monopoles are too close to be seen as separate lumps. For 
SU(2) the constituents are simply found where P{x) = ±12 j38j . 



Fig. 6. The action density for SU(2) calorons with p = f3 and (from left to right) 
M2 = (i.e. trivial holonomy), 0.125 and 0.25 (i.e. TrT'oo = 0), at < = in a plane 
through the constituent locations. 

Important is that when coinciding eigenvalues occur at infinity, that is in 
Voo, some constituent masses will vanish. For trivial holonomy only the nth 
constituent is massive and the action density shows a single lump, the usual 
(deformed) instanton, cmp. Fig. El In the presence of massless constituents 
one cannot take the far field limit, due to the infinite range of the fields in 
these massless cores. One may, however, move massless constituents off to 
infinity. For SU(2) this requires one to take p ^ oo, and one is left with 
a BPS monopole in a singular gauge to adjust the mismatch in boundary 
conditions between caloron and monopole, or put differently, to compensate 
for the magnetic charge that is pushed past infinity 9 . 

As we have explained in the introduction, chiral zero-modes play an ex- 
tremely important role in the physics of instantons. They will also turn 
out to be very useful as a diagnostic tool for studying the properties of the 
caloron solutions. The Atiyah-Index theorem has taught us there is one 
chiral zero-mode for a charge 1 caloron. But when the caloron has "disso- 
ciated" in n constituents the natural question is if the zero-mode is going 
to follow this. There is a compelling reason to expect it will opt for being 



^ In technical terms the locations of constituents are found where log P{x) takes values 
on the boundary of the Weyl chamber |39| . 
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localized on one constituent only. When well separated, the constituents 
become BPS monopoles. These are known to have zero-modes, but only for 
a given range of z values in terms of the Higgs expectation value, as dic- 
tated by the Callias index theorem We will show how this determines 
to which constituent the zero-mode will be localized. Even more useful is 
that we can change this by introducing an arbitrary phase in the boundary 
conditions for the fermions, which reads in the algebraic gauge (Eq. (jSJ) 

+ = exp(27riz/3)Poo^'^(t,x). (51) 

Physical fermions at finite temperature are required to be anti-periodic, 
^'(t, x) = \E'2=i/2/3(i, 2^)- That z determines the location of the zero-mode 
is seen as follows. The gauge transformation g(t) = exp(27rizt) exp(f^g°) 
makes the fermions periodic at the expense of changing = spatial 
infinity to Aq = A'q' — 27riz, cmp. Eq. This acts as an effective mass 
term in the Dirac equation, which differs for each of the gauge components. 
As discussed earlier in this section, each constituent monopole is associated 




Fig. 7. Normalized zero-mode densities at /3 = 1, for the SU{3) caloron of Fig. [S] 
shown for z = jij on equal linear scales (bottom) and on equal logarithmic scales 
for three values of z roughly in the middle of each interval z e (top). 
The zero- mode with anti-periodic boundary conditions is found at z — 30/60. 



with an SU(2) embedding and the two isospin components of these embed- 
dings effectively have "masses" 27r(^m//3 — z) and 27r(^m+i//3 — z). This 
only gives a normalizable zero-mode when these are of opposite sign j4Uj . 
that is for z G [nm/ ^J'm+l/ P]- In the interior of this interval the zero- 
mode is exponentially localized (to the constituent monopole at ym), but at 
z = tJ-m or z = fim+i one of the isospin components becomes massless and 
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the zero-mode delocahzes, having in addition support on the constituent at 
respectively y^-i or Vm+i- 

The expression for the zero-mode density can be given in a simple form 
for arbitrary charge 1 calorons as well |3J. With fJ,m/P^z' <z<fim^i/(3 



fxiz,z') = — — {vm{z')\Ani^i ■ --AiAn ■■■Am- e"^^'*/^|a2i;m(z)) , 

' Vcosh[27r(z-^„/^)r„]y ' ' >^ 4it^ '^•'^^ ' ^' ^ ^ 

where for later use we have introduced also the off-diagonal expression for 
the Green's function fx{z, z') (for z<z' one can use the fact that fx{z', z) = 
f*{z,z')), which can be computed using impurity scattering calculations 
in a piecewise constant potential, see Sec. 12.31 the details of which need 
not concern us here. We illustrate the behavior of the zero-mode in Fig. [7| 
at various values of z for the caloron of Fig. El It is also interesting to 
consider the so-called zero-mode limit (zm) in which x is far removed from 
any constituent location, except the one at ym- This gives^'^ 

pm^ N _ 27r sinh (27rrm(/^m+i//? - z)) sinh (27rrm(z - Hm/P)) 

/3r^sinh(27rr„i.„//3) ' ^^""^ 

At z = \{iJLm + Mm+i)//3 (i.e. z = or z = 1/2/3 for SU(2)) we find 

friz, z) = 7T{rmf3)-^ tanh{7TrmUm/P), (54) 

with — (47r^)~^5^/|™(z, z) giving precisely the zero-mode density of a basic 
BPS monopole, confirming once again the nature of the constituents. In the 
high temperature limit, as long as z 7^ tJ-m/P, the zero-modes become in- 
finitely localized to the constituent locations ym- In this limit the zero- mode 
density is given by (3~^6{x — ym) for fim/P < z < ^m+i/P- The zero-modes 
are therefore ideal for localizing the cores of the constituent monopoles, par- 
ticularly useful for the higher charge calorons discussed below. It should be 
noted that in general the definitions of the constituent locations based on 
the peaks in the energy density, the degeneracy of two eigenvalues in the 
Polyakov loop and the peak in the zero-mode densities will only coincide 
with ym when all constituents have a non-zero mass and are well separated 
(as compared to the size of the monopole cores) from each other. 

Before discussing the higher charge calorons we give the simple ex- 
pression for the SU(2) charge 1 gauge field in the algebraic gauge with 



^'^ Allowing only for errors decaying exponentially in ri^^n, one needs to include 0{r^]^j^) 
dependent shifts in fim, fim+i and Um, see Eq. 1441 and Ref. 31 . 
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Voo = exp(/3a;T3) (hence ^2 = = 

= ^^j^^.'^a^!/ log (/) + |(/>Re ((??^^ - if)ly){Ti + iT2)dyx) , (55) 

where (/•^"'^ = 1 — p-^/^(u;,a;) and x = p'^fx{io, —uj), which has some resem- 
blance to the 't Hooft ansatz in Eq. H13() . It reduces to this ansatz for 
trivial holonomy, w = 0, for which x = ^ ~ 4'~^i one easily checks 
that this gives precisely the Harrington-Shepard solution with (j) = (j)HS: see 
Eq. dJ. In the high temperature limit one finds for non-trivial holonomy 

= 0, (l>^ = {\x-yi\ + \x-y2\ + 1^2 -yi|)/(|^-yi| + \x-y2\ - Im -ml) 

and that = ^fj^^Tsdjy log (f)^ describes a pair of oppositely charged Dirac 
monopoles with the Dirac string on the line connecting them, where (j)^ 
diverges, but outside of which logcf)^ is harmonic. 

3.1. Higher charge calorons 

When ignoring charged components of the gauge field, outside the core 
the Abelian field has unavoidable Dirac strings. We can trace how the 
exact solution instead takes care of the return flux, namely through the 
Abelian component in the magnetic field coming from the commutator of 
the charged components of the non- Abelian field, as is of course well-known 
from the 't Hooft-Polyakov monopole. But from the numerical point of view 




Fig. 8. Approximate superposition of two SU{2) charge 1 calorons (left). The 
logarithm of the action density is plotted through the plane of the constituents 
at t = for (3 — 1 we zoom in by a factor 40 for the transverse direction on the 
would-be Dirac strings (right). 

this will require exponential fine tuning outside the core. It is notoriously 
difficult to find approximate superpositions for magnetic monopoles without 
seeing a remnant of the Dirac string, and in that sense there is "no free 
lunch". Like in the instanton liquid we would like to make approximate 
superpositions of calorons, which also allows us to mix calorons of different 
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charges. However, for those that "dissociate", the constituents should not 
"remember" from which caloron they originated. Although the constituents 
themselves might be well separated, interference can occur in the regions 
between them. An example of this is shown in Fig. |HJ where we added 
two charge 1 caloron gauge fields in the algebraic gauge, called the sum 
ansatz (201 1^ . This preserves the gauge condition but some care is required 
at the gauge singularity, which is removed by a gauge transformation with 
non-trivial winding number (sometimes one refers to these gauge fields as 
being in the singular gauge). Adding the gauge fields after such a gauge 
transformation is performed, destroys the proper decay at infinity. Instead, 
one first smoothly deforms to zero the gauge fields of all other calorons in 
a small neighborhood centered around the gauge singularity one wishes to 
remove. This only costs a small amount of action, particularly when the 
calorons are not too close. In principle a similar construction is possible 
for keeping Dirac strings hidden, but in this case the gauge singularity 
is due to approximating the fields by their Abelian component, far from 
the constituent cores. However, this would require exponential fine tuning. 
Nevertheless, not performing any adjustment the action density along the 
Dirac string (or sheet due to the additional extent in the time direction) 
actually stays finite^^, see Fig. |SJ 

It is due to these complications we felt compelled to analyse the higher 
charge caloron solutions. The disadvantage is that one can only consider 
the exact self-dual solutions this way, and not superpositions of opposite 
charges. On the other hand, it should guarantee absence of Dirac strings, 
as we indeed confirm for a set of exact axially symmetric solutions that share 
all the properties of the charge 1 solutions, as illustrated for the example of 
SU(2) charge 2 solutions in Figs. IHI lllI In the high temperature limit these 
solutions are again described by point-like constituent monopoles . Most 




Fig. 9. Action density on equal logarithmic scales for exact axially symmetric 
SU(2) charge 2 solutions at = 1 for t = 0, on a plane through the constituents. 
Locations can be read off from the dotted lines in Fig. 01 from left to right. For 
the last case see Fig. [TUl 



Only at the gauge singularity the action density of the combined field diverges, which 
may be removed as discussed. 
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of the memory effect has disappeared. In the hght of this it is interesting 
to point out that when further separating the two pairs on the left from 
those on the right, each can be viewed as coming from a single caloron^^. 
However, when the two middle constituents start to get closer than the 
size of their cores, as in Fig. |^ it is more natural to interpret the solution in 
terms of one well "dissociated" caloron formed by the two outer constituents, 
on top of which is superposed a small "non-dissociated" caloron, which is 
no longer static. Indeed, its peak is nearly 0(4) symmetric, whereas the 
other two constituents are static and 0(3) symmetric to a high degree, as is 
appropriate for a BPS monopole (see Fig. |^ (left)). A very subtle memory 




Fig. 10. Action density on equal logarithmic scales for the exact axially symmetric 
SU(2) charge 2 solution with all constituents well-separated (left), the approximate 
solution based on the sum ansatz (middle) and the high temperature limit with 
point-like Dirac monopolcs (right). See also Fig.|51 




Fig. 11. Zero- mode densities on equal logarithmic scales for the caloron in Fig. 1101 
(same horizontal scale). On the left are shown the two periodic zero- modes (z = 0) 
and on the right the two anti-periodic zero- modes {z = 1/2). 

effect, however, remains. It can be shown that for these point-like axially 
symmetric solutions the magnetic charges have to alternate. In particular 
one cannot move one monopole through the other. It may perhaps come as 
a surprise, but in part there is a good reason for being cautious, since as we 
have just seen two oppositely charged constituent monopoles really form a 
small caloron. Since the distance between constituents is given by irp'^/f], 



In this context the parameters and pa appearing in Eq. 14511 can be interpreted as 
the center of mass and size of those calorons. 
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this means that we have to go through a singular caloron when interchanging 
constituent locations on the line. A singular caloron lies on the boundary of 
the moduli (i.e. parameter) space and one cannot use continuity arguments. 

The reason to expect that in general one no longer deals exclusively 
with point-like monopoles comes from the known multi-monopole solutions. 
It is well-known that when putting monopoles of equal charge on top of 
each other, these are deformed in for example the shape of a doughnut j42j . 
as famously illustrated in the scattering of two magnetic monopoles j43j . 
At the technical level this is related to the fact that the Green's function 
mentioned above no longer involves a piecewise constant potential. Only 
for the axially symmetric solutions this was still the case. Nevertheless, for 
SU(2) and charge 2 we have been able to find an exact expression for the 
zero-mode density in the high temperature limit. As long as z 7^ nj, the 
fermions have an infinite mass in this limit, and the zero-modes will vanish 
outside the cores of the constituent monopoles, which themselves need not 
necessarily be isolated points. We thus are able to trace with the help of 
the zero- modes to which region the cores have to be localized. Here we 
will just present the result and discuss its physical significance. For 
lJLrn<z<iJim+i and in the far field limit = —^"^dfVmix), with 

^rn{x) = 7^3^ + — / drdip ^ — , (56) 



27r|f| 47r2 Jr<v \/V 



r 



where y{^) = (Vl — cos ^p, 0, sin 93) , up to an arbitrary coordinate shift 
and rotation. Here P is a scale and k a shape parameter to characterize 
arbitrary SU(2) charge 2 solutions. In this representation it is clear that 
Vm{x) is harmonic everywhere except on a disk bounded by an ellipse with 
minor axes 2PVl-k2 and major axes 2D. Although not directly obvious, 
when k ^ 1 the support of the singularity structure is on two points only, 
separated by a distance 2V. Taking an arbitrary test function /(if) one can 
prove that [HJ 



hm- / /(f)92v„(x)d3x = /(0,0,P)+/(0,0,-P). (57) 

For the caloron k and T) are in general not independent, monopoles 
of different charges have to adjust to each other to form an exact caloron 
solution. As an example we illustrate in Fig. the relation for the case 
studied in Ref. ^J, which is a two parameter family of exact SU(2) charge 
2 calorons solutions in terms of an instanton scale p and relative gauge ori- 
entation angle a. For fixed p it interpolates between two axially symmetric 
solutions. On the left is shown the relation between k and p for some values 
of a. Together with the fact that V = nyflp^iX ^ 0(1 - k^)), this shows 
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that k — > 1 for increasing T>. It can be shown this approach is exponential, 
cmp. Fig. Ha (left). 

The important conclusion is that, when well separated, the constituent 
monopoles become point-like objects. This is a necessary requirement for 




Fig. 12. On the left we plot k versus p for a = (k = 1), a = — tt/IOO and 
a — — 7r/2. On the right are shown the locations of the constituent monopoles (fat 
vs. thin curves for opposite charges) at p = 1/4, varying a from — tt (indicated by 
the arrows) to 0. 

the constituent monopoles to be used as entities to describe the field con- 
figurations at larger distances. Much work of course remains to be done, 
but the results so far have been encouraging. In addition important lattice 
evidence has been accumulated by now, which we will briefiy summarize in 
the next section. 

4. Lattice evidence 

There are a number of lattice studies, that clearly point to constituent 
monopoles to have a dynamical role to play. Two methods have been em- 
ployed. One is based on so-called cooling, the other one uses fermion zero- 
modes. In both cases the purpose is to filter out ultraviolet noise, i.e. one 
is interested in the long distance fiuctuations. 

4.I. Cooling 

Cooling is the process by which one lowers the lattice action through lo- 
cal updates, replacing a link by a suitable combination of neighboring links, 
such that when remaining unchanged it satisfies the lattice equations of mo- 
tion. There are by now many variants, using improved lattice actions (to 
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reduce discretization errors), and criteria to stop the cooling 44j. Such a cri- 
terion is important, because when cooUng too long either all the non-trivial 
fields are removed, or at best one reaches a self-dual solution. The latter is of 
course sometimes done on purpose, so as to reproduce the classical solutions 
on the lattice to quite some precision, or to look for solutions not known 
exactly. For the calorons this has been extensively studied jSH]) even in con- 
nection with a numerical implementation of the Nahm transformation , 
mainly in the presence of so-called twisted boundary conditions |i46) to coach 
the system in having non-trivial holonomy. 

Prom the dynamical point of view, one would like to find how often, and 
with which properties, calorons appear in the long distance fluctuations. In 
this case one would like to start from Monte Carlo generated configurations 
and stop the cooling process when reaching a plateau in the plot for the 
action as a function of the number of cooling steps. The plateau would 
be infinitely long if the configuration is a solution of the lattice equations 
of motion, but would still be sizeable if these are satisfied approximately. 
Therefore, configurations consisting of any number of well-separated instan- 
tons (Q) and anti-instantons (Q) can form such plateaus, with an action 
approximately equal to 87r^(Q -|- Q), whereas the topological charge of the 
equivalent continuum configuration equals Q — Q. An instanton can shrink 
to such a small size p (from the caloron point of view its two constituents 
monolopes getting too close together) that the lattice no longer supports it 
as a solution^^. Also, when instantons and anti-instantons come close to- 
gether they annihilate, this effect is independent of the lattice discretization. 
In both cases the plateau ends, and the cooling curve typically settles down 
at the next plateau where Q + Q has decreased by 1 or 2 unit respectively 
(whereas Q — Q changes by ±1 or respectively). The expectation is that 
the annihilation of oppositely charged configurations is slow, at least when 
they are sufficiently separated, and that from the plateaus one can recover 
the statistical properties of instantons and calorons. 

To prompt the system to have a given holonomy, one can freeze the time- 
like links at the spatial boundary of the box to the required holonomy j48j . 
but in larger volumes it is expected that the confining environment itself 
will provide local regions with a sufficiently coherent background Aq field 
associated with non-trivial holonomy. One important finding has been that 
after cooling the non-trivial holonomy is preserved to some degree I49J . This 
has been analysed for SU(2) in terms of the so-called "asymptotic holonomy" 
Loo, defined as the average of ^TrP{x) over all points x for which the action 
density, summed over t, is smaller than .0001. Histograms of Lqo in the 
confined phase are shown in Fig. (taken from Ref. 011 )• Early on in 



This can be prevented by a suitable choice of improved lattice action |47|. 
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the cooling process a clear peak at Loo = is observed, which becomes a 
flat distribution when cooled down to the plateau associated to Q + Q = 1. 
Important is in particular that the distribution is not becoming peaked 
towards Loo = ±1- Many configurations with well "dissociated" calorons in 
this dynamical setting have been found |49| , and with this method one can 
study the configurations in great detail. For SU(2) one can thus look for the 
monopole centers by finding where the Polyakov loop equals ±12 and use 
the fermion zero-modes for periodic and anti-periodic boundary conditions 
to test if they are localized on the appropriate constituents^^. 





Fig. 13. The "asymptotic holonomy" Loo for an ensemble of 0(200) configurations 
generated on a 16'^ x 4 lattice with periodic boundary conditions for a lattice 
coupling 4/e^ = 2.2. Left is shown the measurement on the first plateau, right 
on the last plateau. For comparison the Haar measure \/l — L^ is shown by the 
dashed curve. Figures taken from Rcf. |49| . 

In these studies one also finds configurations that cannot be directly 
interpreted in terms of instantons, calorons and their possible constituent 
monopoles 1491. Interestingly, there are cases where two constituents appear, 
but of opposite fractional topological charge. These could arise from the 
"annihilation" of two other constituents of opposite duality [SU] but where 
originally each of these, together with one of the surviving constituents, 
formed an (anti-)caloron. More recently these authors also performed cool- 
ing studies for SU(2) at considerably lower temperature than just around the 
deconfining phase transition [HJ. The calorons in this case do not tend to 
"dissociate" in isolated lumps of action density. Nevertheless, constituents 
could be identified through the behavior of the Polyakov loop and were 
characteristic for non-trivial holonomy. This may explain why in the past 
constituent monopoles were never noticed in cooling studies. Finally let us 
mention that many results for SU(3) have now been obtained as well j51j . 

For SU(2) one profits from the fact that z = always lies in the middle of the interval 
[/ii,/i2], whereas z — ^ always lies in the middle of [/i2, 1 + /ii]. 
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4-2. Zero-mode filter 

The use of zero-modes as a filter relies on two observations. The first 
is, as we have seen, that zero-modes quite accurately trace the underlying 
gauge field of instantons, calorons and constituent monopoles. Secondly, a 
zero-mode probes the long distance features of the configuration and ignores 
the ultraviolet high momentum components, otherwise one could never have 
a zero-mode. In some sense the Dirac operator is a particular projection of 
the covariant momentum. For SU(3), comparing periodic and anti-periodic 
fermion boundary conditions 52 , or cycling through all possible phases j53j . 
Eq. (jSlj) . a significantly difi"erent behavior in the two phases was found. In 
the deconfined phase, the proper z-dependent behavior of the trivial holon- 
omy configuration is seen. These are of course the old Harrington-Shepard 
solutions, but their zero-modes were previously only considered for anti- 
periodic boundary conditions. In the confined phase indications for a three 
lump structure is seen (see Fig. 1141 taken from Ref. [SSI) in a reasonable 
fraction of the configurations. These results are based on dynamical con- 





Fig. 14. The zero-mode density for a particular configuration in the confined phase 
on a 6x 20'^ lattice with (Liischer-Weisz) coupling 6/e^ — 8.20, taken from Ref. |53|. 
Shown are slices in {x,y) at {t,z) = (5,9) (left), {t,z) = (2,19) (center) and 
{t, z) = (5, 18) (right), for z = 0.05, 0.3 and 0.65 (from top to bottom). 
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figurations generated with the Liischer-Weisz j54j improved lattice action. 
Configurations that had exactly one chiral zero-mode (which requires the 
use of a chirally improved lattice Dirac operator were singled out. By 
the index theorem, this implies the topological charge is equal to one. The 
dynamically generated configurations can, and typically will, consist of Q+1 
instantons and Q anti-instantons, giving rise to near zero-modes with less 
perfect chiral behavior. Work is in progress to analyse the near zero-modes, 
as well as performing cooling on the SU(3) gauge field configurations [S^ . 
This will help to further rule out the unlikely possibility for the zero-mode 
to jump between instanton, rather than monopole constituents. A more 
theoretical study "ST" has recently shown that it seems indeed impossible 
for the full signature of a single caloron with non-trivial holonomy to be 
emulated by the effect of hopping between instantons. 

Even more remarkable is that a similar structure of multiple lumps for 
the topological charge 1 sector was found on a symmetric lattice of size IG"^ 
and 12"^ at low temperature, well in the confined phase |S31inH]- In this case 
the lattice does not single out one direction as being (imaginary) time. Yet, 
using one of these directions to impose the z dependent boundary condition 
for the fermions, in almost half of the configurations one finds more than one 
lump. The results seem to indicate these lumps are randomly distributed 
over the volume. How to reconcile this with the results found with cooling, 
that pointed to "non-dissociated" calorons at low temperature, is at the mo- 
ment not clear. But let us recall that fractionally charged instantons on the 
torus were long ago introduced by 't Hooft [IHI) and have been extensively 
studied on the lattice [BUI. Twisted boundary conditions and, related to it, 
the quantization of the topological charge in units of 1/n, make it more diffi- 
cult to embed these solutions in a dynamical environment. Nevertheless, on 
the basis of some simple assumptions concerning the dynamical properties 
of these configurations, quasi realistic results have been obtained 61_. It is 
indeed compelling to interpret the 4fcn dimensions of the SU(n) charge k 
instanton parameter space on the torus in terms of the space-time locations 
of kn instanton quarks, as they were called long ago [59 (even though their 
meaning at that time was more abstract). Results for instantons on r^xi?^ 
that in some sense interpolate between the case of calorons and instantons 
on the torus, give further evidence |62j for this conjecture. 



5. Conclusion 

It is clear we can do no justice here to all the results that have been 
obtained in using lattice simulations. The great advantage of using zero- 
modes as a filter is that it minimizes the bias, since it directly uses the 
Monte Carlo configurations themselves. On the other hand, in the cooling 
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studies one has access to the relation between the behavior of the Polyakov 
loop, which plays the role of an order parameter in SU(n) gauge theories, 
and the presence of constituent monopoles, which may give some insight 
in the underlying dynamics. Of course, we would like to make constituent 
monopoles into a precise tool, e.g. for testing the celebrated idea of a dual 
superconductor to describe confinement inSHM]) or as we speculated in the 
introduction, to describe a dense phase of monopoles to lead to deconfine- 
ment of magnetic charges, like quark deconfiment at high density. 

On the more theoretical side we have also not touched upon the role 
these constituent monopoles play in supersymmetric gauge theories |651 167j 
and how the initial motivation for reconsidering calorons with non-trivial 
holonomy came from the study of D-branes {8,. 35, 66^. Nevertheless, we 
feel the applications to the dynamics of non-Abelian gauge theories is very 
promising. We hope to have given the reader some insight in this matter. 
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